Heat Transfer between Weakly Coupled Systems 
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We study the heat transfer between weakly coupled systems with flat interface. We present simple 
analytical results which can be used to estimate the heat transfer coefficient. As applications we 
consider the heat transfer across solid-solid contacts, and between a membrane (graphene) and a 
solid substrate (amorphous Si02). For the latter system the calculated value of the heat transfer 
coefficient is in good agreement with the value deduced from experimental data. 



1 Introduction 

Almost all surfaces in Nature and Technology have 
roughness on many different length scales pj. When two 
macroscopic solids are brought into contact, even if the 
applied force is very small, e.g., just the weight of the 
upper solid block, the pressure in the asperity contact 
regions can be very high, usually close to the yield stress 
of the (plastically) softer solid. As a result good ther- 
mal contact may occur within each microscopic contact 
region, but owing to the small area of real contact the 
(macroscopic) heat transfer coefficient may still be small. 
In fact, recent studies have shown that in the case of sur- 
faces with roughness on many different length scales, the 
heat transfer is independent of the area of real contact Q. 
We emphasize that this remarkable and counter-intuitive 
result is only valid when roughness occur over several 
decades in length scale. 

For nanoscale systems the situation may be very dif- 
ferent. Often the surfaces are very smooth with typically 
nanometer (or less) roughness on micrometer-sized sur- 
face areas, and because of adhesion the solids often make 
contact over a large fraction of the nominal contact area. 
The heat transfer between solids in perfect contact is usu- 
ally calculated using the so called acoustic and diffusive 
mismatch models Q , where it is assumed that all phonons 
scatter elastically at the interface between two materials. 
In these models there is no direct reference to the nature 
of the solid-solid interaction across the interface, and the 
models cannot describe the heat flow between weakly in- 
teracting solids. 

Here we will discuss the heat transfer across perfectly 
flat interfaces. The theory we present is general, but we 
will mainly focus on the case when the interaction be- 
tween the solids is very weak, e.g., of the Van der Waals 
type, as for graphene or carbon nanotubes on many sub- 
strates. We present simple analytical results which can 
be used to estimate the heat transfer coefficient. We 
study heat transfer for solid-solid and solid-membrane 
contacts. We consider in detail the heat transfer be- 
tween graphene and amorphous Si02- For this system 
the calculated value of the heat transfer coefficient is in 
good agreement with the value deduced from experimen- 



FIG. 1: Two solids and 1 in contact. The interfacial surface 
separation is the sum of the equilibrium separation d cq and 
the difference in the surface displacements uo — ui, due to 
thermal movements, where both uo and ui are positive when 
the displacement point along the z-axis towards the interior of 
solid 1. Due to interaction between the solids a perpendicular 
stress (or pressure) ±K(uo(x,t) — ui(x,t)) will act on the 
(interfacial) surfaces of the solid. 



tal data. 
2 Theory 

Consider the interface between two solids, and assume 
that local thermal equilibrium occurs everywhere except 
at the interface. The energy flow (per unit area) through 
the interface is given by[2| 

J = a(T - Ti), 

where To and T\ are the local temperatures at the inter- 
face in solid and 1, respectively. The stress or pressure 
acting on the surface of solid 1 from solid can be writ- 
ten as 

cr(x,t) = K[u (x,t) - tti(x,f)], 

where Uo and u\ are the (perpendicular) surface displace- 
ment of solid and 1 (see Fig. [1]) , respectively, and where 
K is a spring constant per unit area characterizing the 
interaction between the two solids. For weakly interact- 
ing solids the parallel interfacial spring constant Kii is 



usually much smaller than the perpendicular spring con- 
stant K± = K, and we will neglect the heat transfer 
resulting from the tangential interfacial stress associated 
with thermal vibrations (phonons). 
If we define 
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we get 

(7(q,w) = K[u {q,u) - ui(q,w)]. (1) 

Within linear elasticity theory [BJ 

ui(q,w) = Mi(q,u)a(q,u), (2) 

where Mi(q, uj) is determined by the elastic properties 
of solid 1. We consider the heat transfer from solid 
to solid 1. The displacement of an atom in solid is 
the sum of a contribution derived from the applied stress 
— a, and a stochastic fluctuating contribution uof due to 
the thermal movement of the atoms in the solid in the 
absence of interaction between the solids: 

u (q,w) = uof(q,w) - M (q, w)cr(q, w), (3) 

Combining (l)-(3) gives 



ui(q, w) 
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l + X"[Af (q,a;)+Mi(q,w); 



■u f(q,w), (4) 



■Uof(q,w). (5) 



The energy transferred to solid 1 from solid during the 
time period t can be written as 

AE = J d 2 xdt ui(x,t)a(x,t), 

where u = du/dt. One can also write 

AE = (2tt) 3 J d 2 qduj (— iw)« 1 (q, w)<r(-q, -w) 

Using (1), (4) and (5) we obtain 

AE = (2vr) 3 y <2 2 <?doJ 



LoK 2 lmMi(q, u) 



|l + iif[Afo(q,w)+Mi(q,cj)]| 2 



Kf(q,^)| 2 ), (6) 



where we have performed an ensemble (or thermal) av- 
erage denoted by (..). Next, note that (see Appendix 
A) 
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where Aq is the surface area, and 
1 
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is the displacement correlation function. Using the 
fluctuation-dissipation theorem^] we have (see also Ap- 
pendix B and C) 

Cuu (q, w) = ^3 ^ImJ/o (q, w) (8) 

where II (a;) = 7kt> [exp(?kJ / HbTq) — 1] ■ Substituting 
(7) in (6) and using (8) gives the heat current Jo^i = 
AE/A t from solid to solid 1: 

4 r r°° 
Jq^i = j^jz J d q J duj n(w) 



MM (q, oj)lmKMi (q, w) 



|l + if[M (q,^) + Af 1 (q, W )]| 2 ' 

A similar equation with To replaced by T± gives the en- 
ergy transfer from solid 1 to solid 0, and the net en- 
ergy flow J = Jq^,i — Ji^o- The heat transfer coeffi- 
cient a — (Jo->i — Ji^q)/(Tq — T\) gives in the limit 
(T - Ti) 0: 
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(9) 



To proceed we need expressions for A/ (q, w) and 
Mi(q, w). Here we give the M-function for (a) solids, 
(b) liquids and (c) membranes. 

(a) Solids 

For an elastic solid we havejH Q 



M 



i p L (g,^) / w_ 
pc\ S(q,u) \c T 
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where 
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where cl, ct and p are the longitudinal and transverse 
sound velocities, and the mass density, respectively. 

(b) Liquids 

This case can be obtained directly from the solid case 
by letting ct — > 0: 
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(11) 



(c) Membranes 

We assume that the out-of-plane displacement it(x, t) 
satisfies 



Pa 



dt 2 



-kV 2 V 2 u + (t, 



(12) 



where po = fio m o is the mass density per unit area of 
the 2D-system (mo is the atom mass and no the number 
of atoms per unit area), k is the bending elasticity (for 
graphene, k » 1 eV [7j), and cr(x, t) an external stress 
acting perpendicular to the membrane (or rry-plane). Us- 
ing the definition M(q, lu) = u(q,o;)/o'(q, u>) from (12) 
we get 
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Pquj 2 — i0+ 



(13) 



3 Some limiting cases 

Assuming weak coupling between the solids (i.e., K is 
small) and high enough temperature, (9) reduces to 

a = ^3 J d 2 qj^ duj ^Ml m M 1 (q,w)ImA// (q,w). 

""' ' (14) 
In the opposite limit of strong coupling (K — t oo) we get 



„ , °° , dHjuj) ImAf (q,^)ImM 1 (q,^) 
" ■ ' ,UJ dT |Af (q,w) + A/ 1 (q,w)| 2 



(27T) 







(15) 

which does not depend on if. Note also that for very low 
temperatures only very low frequency phonons will be 
thermally excited. Assuming a semi-infinite solid, asu~ 
q — > 0, from (10) we have \M\ ss l/(pcw) — > oo (where 
c is the sound velocity and p the mass density). Thus, 
at low enough temperature (9) reduces to (15) i.e., for 
very low temperatures the heat transfer is independent of 
the strength of the interaction across the interface. The 
physical reason for this is that at very low temperature 
the wavelength of the phonons becomes very long and the 
intcrfacial interaction becomes irrelevant. The transition 
between the two regions of behavior occurs when K\M\k, 
1. Since \M\ k, l/(pcuj) we get K s» pcu. But fiuj ks k B T 
and defining the thermal length At = c/u) = ch/k B T we 
get the condition K w pc 2 / At- Since the elastic modulus 
E w pc 2 we get K « E/Xt. We can define a spring 



constant between the atoms in the solid via k' = Ea, 
where a is the lattice constant. Since K — k/a 2 we get 
k w {a/\^)k' as the condition for the transition between 
the two different regimes in the heat transfer behavior. 
For most solids at room temperature At ~ a, but at very 
low temperatures At >> a which means that even a very 
weak (soft) interface (for which k is small), will appear 
as very strong (stiff) with respect to the heat transfer at 
low temperatures. 

Let us consider the case where the two solids are iden- 
tical, and assume strong coupling where (15) holds. For 
this case we do not expect that the interface will restrict 
the energy flow. If we consider high temperature the ki- 
netic energy per atom in solid will be ~ fee^o so the 
energy density Q « feeTo/cio (where ao is the lattice con- 
stant). Thus if solid 1 is at zero temperature we expect 
the energy flow current across the interface to be of or- 
der J pa Qc/A (the factor of 1/4 results from the fact 
that only half the phonons propagate in the positive z- 
direction and the average velocity of these phonons in 
the z-direction is c/2). Thus we expect a ~ k^cj '(4ag). 
This result follows also from (15) if we notice that for 
Mq = Mi and high temperatures 



(2tt) 3 



d 2 q 
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ImMo(q,u) 
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(16) 



If we assume for simplicity that Mq is given by (11) (but 
the same qualitative result is obtained for solids) then 
the factor involving Mq is equal to unity for lo > c^q and 
zero otherwise. Thus (16) reduces to 



a 
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dq q 
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where we have used that q c « ir/aQ. Thus for identical 
materials and strong coupling (9) reduces to the expected 
result. 

Let us now briefly discuss the temperature dependence 
of the heat transfer coefficient for high and low temper- 
atures. For very low temperatures a is given by (15). 
Consider first a solid in contact with a solid or liquid. 
For these cases it follows that M ~ 1/oj (where we have 
used that u> ~ q) so the temperature dependence of the 
heat transfer coefficient is determined by the term 



Jo 



an(w) 2 



where we also have used that d 2 q ~ oj 2 . Thus we get 



[exp(hu/k B T) - if \k B T 



For high temperatures, and assuming weak coupling, one 
obtains in the same way from (14) that a is tempera- 
ture independent. However, the spring constant (per unit 



area) K may depend on the temperature, e.g., as a result 
of thermally induced rearrangement of the atoms at the 
contacting interface or thermally induced increase in the 
separation of the two surfaces at the interface which may 
be particularly important for weakly interacting systems. 
The temperature dependence of a for the case of solid- 
solid and solid-membrane contacts will be discussed in 
Sec.5. 

4 Phonon heat transfer at disordered interfaces: 
friction model 

At high temperature and for atomically disordered in- 
terfaces, the interfacial atoms will perform very irregular, 
stochastic motion. In this case the heat transfer coeffi- 
cient a can be obtained (approximately) from a classical 
"friction" model. This treatment does not take into ac- 
count in a detailed way the restrictions on the energy 
transfer process by the conservation of parallel momen- 
tum, which arises for periodic (or homogeneous) solids. 
See also Appendix D. 

Let us assume that solid has a lower maximal phonon 
frequency than solid 1. In this case, most elastic waves 
(phonons) in solid can in principle propagate into solid 
1, while the opposite is not true, since a phonon in solid 
1 with higher energy than the maximum phonon-energy 
in solid will, because of energy conservation, be totally 
reflected at the interface between the solids. 

Consider an atom in solid (with mass mo) vibrating 
with the velocity (vm,u_l). The atom will exert a fluctu- 
ating force on solid 1 which will result in elastic waves 
(phonon's) being excited in solid 1. The emitted waves 
give rise to a friction force acting on the atom in solid 
(from solid 1), which we can write as[6j 

F/ = — m 77|| Vj| - m n ± v ± , 

and the power transfer to solid 1 will be 

P = -(Ff ■ v) = m ?7|| <uff > +m i 1± (v'i). 

At high temperatures 

m (v 2 ) = 2fc B T , m (v 2 j_) = k B T . 

Hence 

P = (2 m + V± )k B T . 

A similar formula (with To replaced by T{) gives the 
power transfer from solid 1 to solid 0. Hence 

J = n (2r)\\ + r)±)kB(T - Ti), 

where no = I/ciq is the number of interfacial atoms per 
unit area in solid 0. Thus we get 



For weak interfacial coupling we expect r]± >> m, and 
we can neglect the ny-term in (17). 

The damping or friction coefficient r/± due to phonon 
emission was calculated within elastic continuum me- 
chanics in Ref. |Q|. We have 



■2C> 



VA 



(18) 



where £' w 0.13 (see Ref. [8j). Using that no = l/ a o an< ^ 
substituting (18) in (17) gives 



2 CI 
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where 



-Re 



dx- 
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8tt"\/ " {\-2xf + A{l-x) 1/2 { 1 -x) 1 ' 2 



where 7 = (ct/cl) 2 , and where p — n mo is the (one 
atomic layer) mass per unit surface area of solid 0. There 
are two contributions to the integral One is de- 
rived from the region x < 1 where the integral clearly 
has a non- vanishing real part. This contribution corre- 
spond to excitation of transverse and longitudinal acous- 
tic phonons. The second contribution arises from the 
vicinity of the point (for x > 1) where the denomina- 
tor vanish. This pole contribution correspond to exci- 
tation of surface (Rayleigh) waves. As shown in Ref. 
6], about 65% of the radiated energy is due to the sur- 
face (Rayleigh) phonons, and the rest by bulk acoustic 
phonons. 

We emphasize that (19) is only valid for high tempera- 
tures and weak coupling. A more general equation for the 
heat transfer between solids when the phonon emission 
occur incoherently is derived in Appendix D: 



AA* 
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dn{uj) \mKM Q {u)lmKMi{w) 



dT \l + K[M (u 



Mi (a;)] | 2 ' 
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where 



a = n (2?7|| + r]j_)ki 



(17) 



M(w) = -}- [ d 2 qM(q,tu), 

where the integral is over |q| < q Cl where irq 2 = A*. 
The cut-off wavevector q c is the smallest of qi and q 2 , 
where irq 2 = (2ir) 2 /a^ (where a is the lattice constant) 
and where q 2 = k B T/hcQ (where Cq is the smallest sound 
velocity of solid 0) is the thermal wavevector. For high 
temperatures and weak coupling, for an Einstein model 
of solid 0, (20) reduces to (19) (see Appendix E). 

5 Numerical results 

We now present some numerical results to illustrate 
the theory presented above. We consider (a) solid-solid, 
(b) solid-liquid and (c) solid-membrane systems. 



FIG. 2: The logarithm (with 10 as basis) of the heat trans- 
fer coefficient as a function of the logarithm of the temper- 
ature for weakly interacting solids (soft) with K = 2.52 x 
10 19 N/m 3 , and for solids which interact stronger (stiff) with 
10 times larger K. The dotted line has the slope 3 corre- 
sponding the a ~ T A temperature dependence. 



FIG. 3: The logarithm (with 10 as basis) of the heat trans- 
fer coefficient as a function of the logarithm of the temper- 
ature for weakly interacting solids (soft) with K = 2.52 x 
10 19 N/m 3 , and for solids which interact stronger (stiff) with 
10 times larger K. The solid lines are for coherent phonon 
transmission (from Fig. [2]), and the dashed lines for incoher- 
ent phonon transmission. 



(a) Solid-Solid 

We consider the heat transfer between two solids with 
perfectly fiat contacting surfaces. We take the sound 
velocities and the mass density and the (average) lat- 
tice constant to be that of SiC>2. We consider two 
cases: weakly interacting solids (soft interface) with 
K = 2.52 x 10 19 N/m 3 (see Fig. [5J, and solids with 
stronger interaction (stiff interface), with 10 times larger 
K. In Fig. [2] we show the heat transfer coefficient as a 
function of temperature. Note that for high temperatures 
a is nearly 100 times larger for the stiff case as compared 
to the soft case. This result is expected based on Eq. 
(9b) which shows that a ~ K 2 as long as K is not too 
large or the temperature too low. For low temperatures 
both cases gives very similar results, and for T < 3 K 
the heat transfer coefficient a ~ T 3 . The reason for why 
at low temperatures the heat transfer is independent of 
the strength of the interfacial interaction was explained 
in Sec. 3, and is due to the long wavelength of the ther- 
mally excited phonons at low temperature. 

For incoherent phonon transmission, using (20) we ob- 
tain the result shown by dashed curves in Fig. [3] For 
both the soft and stiff interface the results obtained as- 
suming coherent and incoherent phonon transmission are 
similar. 

(b) Solid-Liquid 

Heat transfer between liquid 4 He and solids was stud- 
ied by KapitzaQ ~ 60 years ago, and R = 1/a is usu- 
ally denoted as the Kapitza resistance 0, [loj]. Let us ap- 
ply the theory to the heat transfer between liquid 4 He 
and a-SiC>2. In the calculation we use a He-substrate 
potential with well depth 10.2 meV and 4 He-substrate 
equilibrium bond distance d cq = 2.2 A which agree with 
the model parameters used in Ref. [l6|. With these 
parameters we get the perpendicular He-substrate vibra- 
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FIG. 4: The calculated contact resistance (multiplied by T 3 ) 
between liquid 4 He and a-Si02 as a function of the tempera- 
ture, with K = 1.18 x 10 19 N/m 3 . 



tion frequency wj_ ~ 91 cm 1 and the spring constant 
K = 1.18 x 10 19 N/m 3 . In Fig. Hwe show the calculated 
contact resistance R (multiplied by T 3 ), as a function 
of the temperature T. In this calculation we have as- 
sumed that all the parameters (e.g., 4 He sound velocity c 
and mass density pn) characterizing the system are tem- 
perature independent [l3T] . The Kapitza resistance has 
been measured (for T > 1 K) for liquid 4 He in contact 
with Quartz [ijj and Sapphire [ll[ and scales roughly with 
temperature as T~ 3 , and the magnitude for T = 1 K is 
roughly 10 times smaller than our calculated result as- 
suming incoherent phonon transfer. Experiment show 
that for T < 0.5 K the Kapitza resistance increase much 
faster with decreasing temperature than expected from 
the R ~ T~ 3 -dependence predicted by our theory and 
most other theories. It is not clear what may be the ori- 
gin of this discrepancy, but it has been suggested to be 
associated with surface roughness. Unfortunately, most 
measurements of the Kapitza resistance was performed 



before recent advances in Surface Science, and many of 
the studied systems are likely to have oxide and unknown 
contamination layers, which may explain the large fluc- 
tuations in the measured contact resistance for nominally 
identical systems. 

(c) Solid-Membrane 

From (13) we get: 



ImMo(q,w) = 7r<5(/i<7 4 — poto 2 ) 



5(u)-ui) (21) 



where uii — (k/po) 1 ^ 2 ? 2 = c (<z)9j where we have defined 
the velocity c(q) = {n/p^Y^q. Substituting (21) in (9) 
and assuming high temperatures so that II (u>) « k^To 
gives: 



Jo- 



2irp, 
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The heat transfer coefficient a = (Jo^i — </i-vo)/(7o — T\) 
is given by 



b f , q „i 
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Using the expression for Mi(q, ui) derived in jH, H and 
uji = c(q)q gives 



where 
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where cl, ct and pi are the longitudinal and transverse 
sound velocities, and the mass density, respectively, of 
solid 1. The cut off wavevector q c sa ir/ai (ai is the lat- 
tice constant, or the average distance between two nearby 
atoms) of solid 1 . 

There are two contributions to the integral £. One is 
derived from c(q) > cl, but for graphene on a-SiC>2, this 
gives only ~ 10% of the contribution to the integral. For 
c{q) < cl the term after the Re operator is purely imag- 
inary (and will therefore not contribute to the integral), 
except for the case where the denominator vanish. It 
is found that this pole-contribution gives the main con- 
tribution (~ 90%) to the integral, and corresponds to 
the excitation of a Rayleigh surface (acoustic) phonon of 
solid 1. This process involves energy exchange between a 
bending vibrational mode of the graphene and a Rayleigh 



surface phonon mode of solid 1. The denominator vanish 
when c(q) — cr where 
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= 0. 



Note that the Rayleigh velocity cr < ct but close to ct- 
For example, when cl/ct = 2, cr w 0.93ct, and the pole 
contribution to the integral in £ is 0.083. Note that (23) 
is of the same form as (19), and since £' ss 0.13 f=s ^ they 
give very similar results. 

In the model above the heat transfer between the solids 
involves a single bending mode of the membrane or 2D- 
system. In reality there will always be some roughness at 
the interface which will blurred the wavevector conserva- 
tion rule. We therefore expect a narrow band of bending 
modes to be involved in the energy transfer, rather than 
a single mode. Nevertheless, the model study above as- 
sumes implicitly that, due to lattice non-linearity (and 
defects), there exist phonon scattering processes which 
rapidly transfer energy to the bending mode involved in 
the heat exchange with the substrate. This requires very 
weak coupling to the substrate, so that the energy trans- 
fer to the substrate is so slow that the bending mode 
can be re-populated by phonon scattering processes in 
the 2D-system, e.g., from the in-plane phonon modes, 
in such a way that its population is always close to what 
would be the case if complete thermal equilibrium occurs 
in the 2D-system. This may require high temperature in 
order for multi-phonon scattering processes to occur with 
enough rates. 

We now consider graphene on amorphous Si02- 
Graphene, the recently isolated 2D-carbon material with 
unique properties due to its linear electronic dispersion, 
is being actively explored for electronic applications 14 1 . 
Important properties are the high mobilities reported es- 
pecially in suspended graphene, the fact that graphene is 
the ultimately thin material, the stability of the carbon- 
carbon bond in graphene, the ability to induce a band 
gap by electron confinement in graphene nanoribbons, 
and its planar nature, which allows established pattering 
and etching techniques to be applied. Recently it has 
been found that the heat generation in graphene field- 
effect transistors can result in high temperature and de- 
vice failure 15]. Thus, it is important to understand the 
the mechanisms which influence the heat flow. 

The graphene — a-Si02 interaction is probably of the 



Van der Waals type. In Ref. [17( the interaction between 



the graphene C-atoms and the substrate Si and O atoms 
was assumed to be described by Lennard- Jones (LJ) pair- 
potentials. Here we use a simplified picture where the 
substrate atoms form a simple cubic lattice with the lat- 
tice constant determined by a\ — (m/pi) 1 / 3 0.25 nm, 
where m = (msi + 2m )/3 « 3.32 x 10 -26 kg is the 
average substrate atomic mass, and p\ w 2200 kg/m 3 
the mass density of a-SiC>2. We also use the effective 
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FIG. 5: The calculated graphene — a-Si02 interaction energy 
17(d) per graphene carbon atom, as a function of the separa- 
tion d (in nm) between the center of a graphene carbon atom 
and the center of the first layer of substrate atoms. See text 
for details. 
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FIG. 6: The logarithm of the contact resistance for graphene 
on a-Si02 as a function of the logarithm of the temperature. 
Square symbols: measured data from Ref. [2(j. Solid lines: 
the calculated contact resistance using Eq. (24) with hcoo = 
11 meV (upper curve) and Hloq = 7 meV (lower curve). 



LJ energy parameter, e = (egj + 2eo)/3 ~ 5.3 meV, and 
the bond-length parameter a = (<Tsi + 2cro)/3 » 0.31 nm. 
With these parameters we can calculate the graphene— a- 
Si02 interaction energy, U(d), per graphene carbon 
atom, as a function of the separation d between the center 
of a graphene carbon atom and the center of the first layer 
of substrate atoms. We find (see Fig. [SJ the graphene— a- 
Si02 binding energy = —U(d e q) ~ 35 meV per car- 
bon atom, and the force constant k = U"(d cq ) (where 
d cq ~ 0.32 nm is the equilibrium separation) k = Ka^ = 
2.4 N/m per carbon atom. This gives the perpendic- 
ular graphene — a-Si02 (uniform) vibration frequency 
lux_ sd (fc/mn) 1 / 2 ~ 55 cm -1 , which is similar to what is 
observed for the perpendicular vibrations of linear alkanc 
molecules on many surfaces (e.g., about 50 — 60 cm -1 for 
alkanes on metals and on hydrogen terminated diamond 
C(lll)l8]). Using K = kja\ = 1.82 x 10 20 N/m 3 , and 
the transverse and longitudinal sound velocities of solid 

1 (ct = 3743 m/s and cl = 5953 m/s), from (19) we 
obtain a » 3 x 10 8 W/Km 2 . 

The heat transfer coefficient between graphene and a 
perfectly flat a-Si02 substrate has not been measured di- 
rectly, but measurements of the heat transfer between 
carbon nanotubes and sapphire by Maune et al[19j in- 
dicate that it may be of order a w 8 x 10 8 W/m 2 K. 
This value was deduced indirectly by measuring the 
breakdown voltage of carbon nanotubes, which could be 
related to the temperature increase in the nanotubes. 
Molecular dynamics calculations 17] for nanotubes on 
a-Si02 gives a « 3 x 10 8 W/m 2 K (here it has been 
assumed that the contact width between the nanotubc 
and the substrate is 1/5 of the diameter of the nan- 
otubeh Finally, using a so called 3 ui method, Chen 
et al[20j have measured the heat transfer coefficient a « 

2 x 10 8 W/m 2 K. 

We now discuss the temperature dependence of the 
heat transfer coefficient. If we assume that most of the 
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FIG. 7: The frequency of the graphene bending mode be- 
comes equal to the frequency of the Rayleigh mode when 
Wbend(<?) = ujR(q) ■ The Rayleigh mode dispersion was mea- 
sured for a-quartz (0001) [241] and the bending mode disper- 
sion was calculated using u;bcnd = (ft/po) 1 ^ 2 *^* with the bend- 
ing stiffness k = 1.1 eV as obtained in Ref. [7]]. 



heat transfer is via a substrate phonon mode at the fre- 
quency wo, then the temperature dependence of a should 
be given by 



dT 



x 2 e x 



(e x - l) 2 



(24) 



where x = TiuiQ/k^T. In Fig. [6] we show the tempera- 
ture dependence of the heat transfer coefficient measured 
by Chen et al [20] for 42 K < T < 310 K. The solid 
lines has been calculated using (24) with Tluiq = 11 meV 
(upper curve) and 7 meV (lower curve). In our model 
all the vibrational modes have linear dispersion relation, 
e.g., huj — cr(7 for the Rayleigh mode, and the frequency 
where the graphene bending mode become equal to the 
frequency of the Rayleigh mode will occur at higher fre- 
quency than expected using the measured Rayleigh mode 
dispersion relation. This is illustrated in Fig. [7] where 



we show the measured Rayleigh mode dispersion for a- 
quartz (0001) [24J. Note that the frequencies of the bend- 
ing mode and the Rayleigh mode become equal when 
Tilo r; 7 meV. However, using this excitation energy in 
(24) gives too weak temperature dependence. There are 
two possible explanations for this: 

(a) In an improved calculation using the measured dis- 
persion relations for the substrate phonon modes, emis- 
sion of bulk phonons may become more important than in 
the present study where we assumed the linear phonon- 
dispersion is valid for all wavevectors q. This would make 
higher excitation energies more important and could lead 
to the effective (or average) excitation energy 11 meV 
necessary to fit the observed temperature dependence. 

(b) As pointed out in Sec. 3, the model developed 
above for the heat transfer involves a single, or a narrow 
band, of bending modes of the membrane or 2D-system. 
In order for this model to be valid, the coupling to the 
substrate must be so weak that the energy transfer to the 
substrate from the bending mode occurs so slowly that 
the mode can be re-populated by phonon scattering pro- 
cesses, in such a way that its population is always close to 
what is expected if full thermal equilibrium would occur 
within the 2D-system. This may require high temper- 
ature in order for multi-phonon scattering processes to 
occur by high enough rates. This may contribute to the 
decrease in the heat transfer coefficient observed for the 
graphene — a-SiC>2 system below room temperature j20| . 

Finally, we note that recently it has been suggested flBL 
I22] ] that the heat transfer between graphene and a-SiC>2 
may involve photon tunneling [2^]. That is, coupling via 
the electromagnetic field between electron-hole pair exci- 
tations in graphene and optical phonons in a-SiC>2. How- 
ever, our calculations indicate that for graphene adsorbed 
on a-SiC>2 the field coup ling gives a negligible contribu- 
tion to the heat transfer 2l| . 



6 Role of surface roughness 

Surface roughness has usually a strong influence on the 
heat transfer between macroscopic solids 0. For hard 
solids the area of real contact may be very small com- 
pared to the nominal contact area, and in these cases 
most of the heat may flow in the air film separating the 
non-contact region. The heat transfer via the area of 
real contact is determined not just by the heat transfer 
resistance across the contacting interface (of atomic scale 
thickness) as studied above, but often most of the heat 
flow resistance is caused by the so called spreading resis- 
tance, related to the interaction between the heat flow fil- 
aments which emerge from the areas of real contact. This 
latter contribution depends on the wide (fractal-like) dis- 
tribution of surface roughness length scales exhibited by 
most surfaces of macroscopic solids, see Fig. [5J One can 
show that the total heat transfer resistance is (approxi- 




FIG. 8: Heat flow in the contact region between a rigid block 
with a flat surface (bottom) and en elastic solid with a ran- 
domly rough surface (top). The orange lines denote the heat 
current flux lines in the upper solid. The heat current fila- 
ments expand laterally until the filaments from the different 
contact regions touch each other. The "interaction" between 
the filaments gives rise to the spreading resistance. Because 
of the fractal nature of most surfaces the interaction between 
the heat flow filaments occur on many different length scales. 



mately) the sum of the two mentioned contributions: 
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where 1 / a S pred is the spreading resistance studied in Ref. 
0, and l/a c the resistance which determines the tem- 
perature jump (on atomistic length scale) across the area 
of real contact. One can show that (see Appendix F): 



1 1 



d 2 x J z 2 (x) 



where J z (x) is the heat current at the interface, Jo the 
average heat current, and ab the (boundary) heat trans- 
fer coefficient studied above (see Sec. 2). If the heat 
current would be constant through the area of real con- 
tact, then J z = (Ao/A)Jq, where A is the area of real 
contact. In this case we get a c ~ {A/A^a^ and 
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a 
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^sprcd 



l_Ao 

ctb A 



(25) 



For most hard macroscopic solids the local pressure in 
the contact regions is very high, which may result in 
"cold welded" contact regions with good thermal con- 
tact, in which case the contribution from the spreading 
resistance dominates the contact resistance. However, for 
weakly coupled microscopic solids the contribution from 
the second term in (25) may be very important. 

In Ref. 15| the temperature profile in graphene under 
current was studied experimentally. The heat transfer co- 
efficient between graphene and the a-Si02 substrate was 



determined by modeling the heat flow using the standard 
heat flow equation with the heat transfer coefficient as 
the only unknown quantity. The authors found that us- 
ing a constant (temperature independent) heat transfer 
coefficient a pa 2.5 x 10 7 W/m 2 K the calculated temper- 
ature profiles in graphene are in good agreement with 
experiment. This a is about 10 times smaller than ex- 



pected for perfectly flat surfaces (see Sec. 5(c)). In 21] 
we have studied the heat transfer between graphene and 
a-Si02- We assumed that because of surface roughness 
the graphene only makes partial contact with the SiC>2 
substrate, which will reduce the heat transfer coefficient 
as compared to the perfect contact case. The analysis 
indicated that the spreading resistance contribution in 
(25) may be very important, and could explain the mag- 
nitude of the observed heat contact resistance. However, 
assuming that (due to the roughness) A/Aq pa 0.1, the 
second term in (25) becomes of the same order of mag- 
nitude as the measured heat contact resistance. Thus, in 
this particular application it is not clear which term in 
(25) dominates the heat resistance, and probably both 
terms are important. 



value deduced from the experimental data. 
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Appendix A 

Here we prove equation (7). We get 



(Kf(q^)| 2 } = 



1 



(2tt) 6 



d 2 xdtd 2 x'dt' 



X (U0i (X, t) U0{ (x', t'^Mx-x')-^-*')] 



7 Summary 

To summarize, we have studied the heat transfer be- 
tween coupled systems with flat interface. We have pre- 
sented simple analytical results which can be used to 
estimate the heat transfer coefficient. The interaction 
between the solids is characterized by a spring constant 
(per unit area) K. The formalism developed is general 
and valid both for strongly interacting (K — > oo) and 
weakly interacting (K — > Q) solids. We have shown that 
at low enough temperatures, even a very weak interfacial 
interaction will appear strong, and the heat transfer is 
then given by the limiting formula obtained as K — > oo. 
Earlier analytical theories of heat transfer Q does not 
account for the strength of the interaction between the 
solids, but correspond to the limiting case K — > oo. How- 
ever, we have shown that at room temperature (or higher 
temperatures) the heat transfer between weakly inter- 
acting solids may be 100 times (or more) slower than 
between strongly interacting solids. 

Detailed results was presented for the heat transfer be- 
tween a membrane (graphene) and a semi-infinite solid 
(a-Si02). For this case the energy transfer is dominated 
by energy exchange between a bending vibrational mode 
of the graphene, and a Rayleigh surface phonon mode of 
the substrate. This model assumes implicitly that, due 
to lattice non-linearity (and defects), there exist phonon 
scattering processes which rapidly transfer energy to the 
bending mode involved in the heat exchange with the 
substrate. This may require high temperature in order 
for multi-phonon scattering processes to occur at high 
enough rate. The calculated value of the heat transfer 
coefficient was found to be in good agreement with the 
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Appendix B 

Here we present an alternative derivation of Eq. (8). 
Assume that the two solids interact weakly. In this case 
the energy transfer from solid to solid 1 is given by (6) 
with K ->• 0: 

AE= (2tt) 3 J d 2 qdoj w^ImM^q, w)([it of (q, w)\ 2 ). 

(Bl) 

At thermal equilibrium this must equal the energy trans- 
fer from solid 1 to solid given by 

AE = (2vr) 3 J d 2 qduj wAT 2 ImM (q, w)(K f (q, w)| 2 ). 



From (Bl) and (B2) we get 



(B2) 



ImM 1 (q,a;)(|wof(q,^)| 2 ) = ImAf (q,o;)(|w lf (q,a;)| 2 ) 

(B3) 



We now assume that solid 1 is a layer of non-interacting 
harmonic oscillators. Thus if pi is the mass per unit area 
we have 

p\U\ + u\u\ = a 



or 



Thus 



/ \ o-(q,w) 

Ml,") = — ,2 _ , , 2 \ 



Pl (uj 2 - UJ 2 ) - l0+ 



Mi(q,w) - 



pi(w? -u; 2 ) -t0+ 



and for cj > 



Im¥i(q,w) = 5(cj - wi) (S4) 

2/9iWi 



We write Mi in the standard form: 



1 



1/2 



x (V i(q - X - Wlt) + 6+e- t(q ' 
so that for uj > 0: 



X — Wit) 



ui(q,w) 



Thus we get 



1 / ft 



(2tt) 2 



1/2 



6 q <5(w — wi) 



where we have used that 



1 f t 

[S(w - uj!)} 2 = S(uj - wi)— / dt = 5(w - wi) — 

Using that 

(& q fe++&+fo q ) - [2n( Wl )+l](2^) 2 (5(q-q) = [2n(wi)+l]A 
we get 

(l-if(q,-)| 2 ) = ff B ^ (»H + ^ % - «i) 

(B5) 

Combining (B3)-(B5) gives 



(l " l,r(q ' J,)|a> = w(" ( '" ) ^ ,I,,1 ' U " (q ^ <L,(il 



Eq. (8) is a standard result but the derivation is re- 
peated here for the readers convenience. Let us write the 
Hamiltonian as 

H = H a + J d 2 x u(x,t)a(x,t) 

where <r(x, t) is an external stress acting on the surface 
z = of the solid. We first derive a formal expression for 
M(q, uj) defined by the linear response formula 

(w(q,cj)) = M (q, w)<r(q, w) 

We write (u) — Ti(pu) where the density operator satis- 
fies 



>(t-t')/h 



We write p = p + p\ and get 

Thus using (u) = Tr(piu) we get 
(u) = ^rj d 2 x'dt' 6(t-t')([u{xL,t),u{x',t')])a(x.',t') 

= hj d2;r ' rfi ' 6 '(*- < ')([w(x-x',t-t'),w(O,0)])cr(x',t') 



Thus 



M(q,w) = 4r / rf 2 a;* 0(t)([u(x, i), tt(0, 0)])e 
in 7 



-{6 q 6+ + 6+& q )(5(cj - wi) whcrc 



-ifq-x— u;i) 

(Ci) 



u(x,t) = e- lHot / ?l u(x,0)e 



iHat/h 



Let |n) be an eigenstate of corresponding to the 
energy E n . Using (CI) we get 



M(q 
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-i(E n -E m )t/h 
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Appendix C 



/(n|u(x,0)|m)(m|u(0,0)|n) 



(n|w(0, 0)|to)(to|w(x, 0)|n) 
+ i(.E ri -.E m )/7i + zw-0+ 



(C2) 



where /3 = l/fc B T and Z = £ n exp(-/3£„). From (C2) 
we get 



ImM(q,w) 



e ^"e- lq ' x 



x «n|u(x, 0)|m><m|«(0, 0)|n)(-J[w - (£?„ - E m )/K\) 

+(n\u(0,0)\m)(m\u{x,0)\n)5[u; + (E n - E m )/h}) (C3) 

Changing summation index from (n, m) — >• (m, n) the 
(to, n)-dependent part of the second term in (C3) can be 
rewritten as 

^e- /3B " l (m|u(O,0)|n)(n|u(x,0)|m) 



xJ[w-(£^-f^)/R] 



-PE np -P{E m -E n ) 



n|u(x, 0)|m)(m|u(0,0)|n) 



^ e- /3 - E "(n|u(x, 0)|m)(m|u(0, 0)|n) 



x6[u-(E n -E m )/h] 

Replacing the second term in (C3) with this expression 
gives 



ImM(q, W ) = i(e^-l) 



x J d 2 xdt e- l ^-^(u(^t)u(O,0)) 



= 2^(^-1) (27r) 3 C m (q,c) 

From the last equation follows the fluctuation-dissipation 
theorem: 



Appendix D 

We assume high temperatures and interfacial disorder. 
In this case the elastic waves generated by the stochastic 



pulsating forces between the atoms at the interface give 
rise to (nearly) incoherent emission of sound waves (or 
phonons). Thus, we can obtain the total energy trans- 
fer by just adding up the contributions from the elastic 
waves emitted from each interfacial atom. Assume for 
simplicity that the interfacial atoms of solid for a sim- 
ple square lattice with lattice constant a$. Consider the 
atom at x = and let uo(t) denote the vertical displace- 
ment of the atom. The force 

F(t) =fc[«o(*)-«i(t)], 



or 



F(uj) = k[u (w) - ui(w)], 



(£>1) 



is acting on solid 1 at x = 0. We can write k — Kclq 
where K is the force constant per unit area (see Sec. 2). 
The force F(t) gives a stress 

cr(x,t) = F(t)5(x) 

acting on solid 1. We can also write 

( r(q, W ) = (27r)- 2 FH 

Note that 

Ui(u}) = iti(x = 0,w) = J d 2 q Ui(q, uj) 
= / d 2 q Mi(q,w)a(q, uj) 



^2 / < /J " ^(q.^FU) 



(D2) 



where 



Mi(w) 



d 2 q Mi(q,w) 



(2n) 2 

In a similar way one get 

uo(w) = uof (w) - M (u)F(u). (D3) 
Combining (D1)-(D3) gives 



, x 1 + fcMiH 

MoM = i + mmoM + ^M] ^- (D5) 

The energy transferred to solid 1 from solid during the 
time period to can be written as 



AE = N 



J dt in(t)F(t), 



where N = Ao ja\ is the number of interfacial atoms of 
solid 0. One can also write 

AE = 2-kN j duj (-iw)ui(u))F(-u) 



Using (Dl) and (D4) and (D5) we obtain 
wPlmMi(id) 



AE = 2nN 



+ fc[M (w) + MiH]| 



2 (KfHI 2 >, 



(D6) 

where we have performed an ensemble (or thermal) av- 
erage denoted by (..). Next, note that 



(KfMI 2 ) 
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dtdt' (u 0f (t)u f(t')}e~ 4aj( *~*' ) 



(2tt) 2 
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where 



^ j dtdt'(u of (t)u of (0))e iut = ^C uu (u;), (D7) 



C uu (uj) i- j (» ()r (/)» () r(().))« 



is the displacement correlation function. Note that 

C uu {uj) = j d 2 q C utt (q,w) 
Thus, using (8) we get 

Cuu(u) = — ^r^^Irn / d 2 q M (q,w) 



2 n(w) 

27T U> 



ImMo(w) 



(D8) 



Substituting (D7) in (D6) and using (D8) gives the heat 
current Jo^i = AE/A to from solid to solid 1: 



Jo-¥i — 



AA* 



a* r 



du H(ui)- 



lmK Mq{oj)1vd.K Mi(u>) 



where A* = (2ir) 2 ja\ is the area of the Brillouin zone 
and where we have defined 

M(w) = a 2 Q M(w) = 1^1 d 2 q M(q, w) 

where the q- integral is over |q| < q c with irq 2 = A*. A 
similar equation with T replaced by T\ gives the energy 
transfer from solid 1 to solid 0, and the net energy flow 
J = Jo^i — </i->o- The heat transfer coefficient a = 
(Jo-n - Ji^o)/(T -Ti) gives in the limit (T -Ti) -> 0: 



4A* r 



dijj 



dn(oj) ImKMoi^lmKMxjuj) 

(D9) 



The derivation above is only valid for high temperature 
where k B T > Tiwo, where fi-wo is the highest phonon en- 
ergy of solid 0. However, we can apply the theory (ap- 
proximately) to all temperatures if we take the cut-off 
wavevector q c to be the smallest of qi and q2, where 
irq 2 = (2ir) 2 /a 2 , (where ao is the lattice constant) and 
where 02 = h-aT/Tico (where Co is the smallest sound ve- 
locity of solid 0) is the thermal wavevector. 

Appendix E 

Here we show that (20) reduces to (19) for high tem- 
peratures and when solid is described by an Einstein 
model. At high temperatures and weak interfacial cou- 
pling, (20) becomes 

Ah-aA* f°° 

a= / dw ImA'Mo(w)ImA'Mi(w) (El) 

We assume for solid that 
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M (u) 



Po(wq -w 2 ) -«0+ 



where p n = mo/a 2 , is the mass per unit area. Thus 



ImMo(w) = 



2p LO 

Substituting this in (El) gives 
'cb 
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ImMi(uo) 



kBK 2 1 / , 2 
(2n) 2 p a ujQ 
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d 2 q ImMi(q, w ) 



Substituting () into this equation and denoting q 
(loq/c^x 1 / 2 gives 

k B K 2 



x —Re 



dx- 



2 (7 - x) 



1/2 



Jo ' (1 - 2x) 2 + 4 (1 - x) 1 ' 2 (7 - x) 
where 7 = (ct/cl) 2 , which agree with (19). 
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(E2) 



Appendix F 

Let l/a c be the interfacial contact resistance associ- 
ated with the jump in the temperature (on an atomistic 
scale) in the contact area between two solids, and let 
1 / a sp rcd be the spreading resistance associated with the 
interaction between the heat filaments emerging from all 
the contact regions. Since these two resistances act in 
series one expect the total contact resistance to be the 
sum of the two contributions, i.e., 



1 

a 



1 1 
1 

^sprcd 



We can prove this equation and Eq. (24) using the 
formalism developed in Ref. Q. We assume that all 
the heat energy flow via the area of real contact. In 
this case the interfacial heat current J z (x) vanish in 
the non-contact area. In the area of real contact the 
temperature T(x, z) change abruptly (on an atomistic 
scale) when z increases from z = — + (in solid 0) to 
z = + (in solid 1), and the jump determines the heat 
current: J z (x) = a h [T(x, -0) - T(x, +0)]. If we denote 
■0(x) = T(x, -0) - T(x, +0) the equation 

J a (x)[J,(x)-ObV(x)] =0 

will be valid everywhere at the interface. From this equa- 
tion we get 

J d 2 q'J z (q - q')[J.(qO - «bV>(q')] = 

Following the derivation in Sec. 2.2.1 in Ref. 0] we get 
instead of Eq. (20) in Ref. Q the equation 

a k A Jfi J q 

+-7^72 / ^ M = 7^— + - 
ab A J J a sp rod a c 

Here Jo is the average or nominal heat current, AJ z (x) = 
J z (x) — Jo, and k an effective heat conductivity (k _1 = 
Kq" 1 + K^f 1 ). The first term in (Fl) is the spreading re- 
sistance studied in Ref. Q while the second term is the 
contribution from the temperature jump on the atomistic 
scale accross the area of real contact. 
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